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1^ ■ The pionic decay rates of the excited L = 0,1 D mesons are calculated with a Hamiltonian 

model within the framework of the covariant Blankenbecler-Sugar equation. The interaction be- 
tween the light quark and charm antiquark is described by a linear scalar confining and a screened 
one-gluon exchange interaction. The decay widths of the D* mesons obtain a contribution from 
^^' ' the exchange current that is associated with the linear scalar confining interaction. If this contri- 

f^ , bution is taken into account along with the single quark approximation, the calculated decay rates 

f^ ' of the charged D* mesons are readily below the current empirical upper limits if the axial coupling 

constant of the light constituent quarks is taken to be g\ — 0.87, but reach the empirical upper 
limits if g\ = 1. With the conventional values for g\, the calculated widths of the Di and D2 
mesons fall somewhat below the experimental lower limits, leaving room for other decay modes as 
f^ {_ well, such as tttt decay. The unrealistically large contribution from the axial charge operator to the 

Qj ■ calculated pion decay width of the Di meson is suppressed by taking into account the exchange 

charge effects that are associated with the scalar linear confining and vector one-gluon exchange 
interactions. The predicted values for the pionic widths of the hitherto undiscovered L = 1 Dl 
and Do mesons are found to be smaller than previous estimates. 
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1 Introduction 

The pion decay rates of the excited charmed mesons - the D mesons - may provide direct information 
on the strength of the pion coupUng to Ught constituent quarks. As the charm quark in the D mesons 
does not couple to pions, the decay mechanism is determined by the pion couphng to the hght flavor 
constituent quark. A first assumption is that this couphng is independent of the interaction between 
the hght quark (or antiquark) and the charmed antiquark (or quark). While this may be considered 
as a satisfactory approximation for the axial current part of the pion-quark coupling, it leads to large 
overestimates of the decay widths in the case of the axial charge term, a problem that may be cured 
by the two-quark mechanism that is associated with the confining interaction between the light and 
charmed quark in the qc (qc) system. 

Because of the large velocities of the confined quarks in D mesons, and in view of the small mass of 
the light constituent quarks, the D mesons have to be described as relativistic interacting two-particle 
systems. The simplest way to achieve such a description is to employ a covariant three dimensional 
reduction of the Bethe-Salpeter equation and a corresponding quasipotential representation of the 
interaction between the constituents. The hitherto considered quasipotential descriptions for the D 
mesons are based on the Gross m and the Blankenbecler-Sugar equations M, m|. Both approaches 
have been shown to yield reasonable predictions for the D meson spectrum with combinations of linear 
confining and one-gluon exchange interactions between the quark and antiquark H, pi pi . 

This dynamical model, with static interactions in the qQ system applied in the solution of the 
Gross equation, has recently been used to calculate the pseudoscalar meson decay rates of the D 
mesons B. A fair description of the decay rates of the D* mesons was achieved, under the assumption 
that the pions couple to the light constituent quarks by the standard (pseudovector) coupling. On the 
other hand, the ratio of the widths of the D2(2460) and the 151(2420) mesons was found to be ~ 2.5, 
while the empirical ratio is ~ 1.3. In view of the importance of determining the form and strength of 
the pion-quark coupling, an analogous calculation within the framework of the Blankenbecler-Sugar 
equation is performed here. There is no obvious reason for preferring one or the other quasipotential 
framework, besides that of calculational convenience. The Gross equation framework has the virtue 
of reducing to a Dirac equation for the light quark in the infinite mass limit of the heavy (charm) 
quark. The Blankenbecler-Sugar equation has an interpretational advantage in its formal similarity to 
the standard Schrodinger equation framework. 

To calculate the pion decay rates of the excited D mesons, the wavefunctions that have been 
obtained in ref. [|6| by solving the Blankenbecler-Sugar equation in configuration space are employed in 
order to achieve a description of the states in the D meson spectrum. These wavefunctions correspond 
to a model which describes the interaction in the qQ system as a scalar linear confining interaction 
combined with a screened relativistic one-gluon exchange interaction, which yields a spectrum that 
agrees with the empirically known part of the spectrum. The model leads to hyperfine splittings that 
agree well with recent NRQCD lattice calculations in the quenched approximation |g[. The rates for the 
decays of the form D' -^ Dtt are then obtained by calculating the matrix elements of the pion creation 
operator between the excited and ground states using such wave functions. Because of the small mass of 
the light flavor quark, which couples to the pions, the non-local structure of the pseudoscalar coupling 
of the pion to the quark has to be treated in unapproximated form. The static (local) approximation 
to this vertex function is shown to imply an overestimate by about a factor 2. 

The present empirical information on the widths of the excited charm mesons remains very incom- 
plete. Absolute values, with large uncertainty ranges, are known for the decay widths of the 151(2420) 
and 1)2(2460) mesons, but for the i:'*(2010)='= and the £)*(2007)° only upper hmits are available at 



present Q . The pion decays of the D* mesons to D mesons are P— wave decays generated by the axial 
current operator. If the pion is assumed to be emitted by a single quark operator with the conventional 
value for the pion-quark coupling constant, the calculated widths for the decays D* -^ Dtt fall well 
below the present empirical upper limits in the case of the charged D* mesons (The empirical upper 
limit on the width of the neutral D* meson is too large to be constraining). Upon addition of the 
contribution from the two-body axial exchange (pair) current that is associated with the linear scalar 
confining interaction, the calculated widths reach the empirical upper limits in case of the charged D* 
mesons, if the value of the axial coupling constant of the light constituent quarks is taken to be g^ = 1. 

The single quark mechanisms for pion production lead to a considerable overprediction of the S'-wave 
pion decay widths of the Di mesons. This overestimate may be reduced by invoking the two-quark 
mechanism, which is naturally associated with the scalar confining and vector one-gluon exchange 
interactions. Consequently, the predicted widths of other D meson states that decay by an S'-wave 
mechanism are also suppressed by large factors. This effect was first hinted at in ref. [g|. In the case of 
the scalar confining interaction the simplest description of this mechanism is to view it as an effective 
increase of the constituent quark mass from rriq to mq + cr, where c is the confining string tension. 
Since the quark mass appears in the denominator of the transition amplitude for S-wave pion decay, 
this increase of the constituent quark mass leads to a large reduction of the associated matrix element. 
An analogous suppression of S'-wave pion decay modes was achieved in ref. [|| as a correction to the 
one-quark operator through coupling to negative energy states. This two-body mechanism is analogous 
to that, which is required for a realistic description of the Ml decays of charmonium and heavy light 
mesons ||, |7| . 

Once this two-quark mechanism is taken into account, the calculated pion decay widths of the 
Di (2420) and Z?2 (2460) mesons fall somewhat below the empirical values if the axial coupling of the 
constituent quark is taken to be less than 1. If, on the other hand, the matrix element of the pion 
decay amplitude is evaluated in the non-relativistic approximation, the calculated decay widths exceed 
the empirical values, in agreement with the result of ref. m. In the case of the Di meson, the S— wave 
pion decay mode is found to contribute significantly, so that in the end the net ratio of the calculated 
widths of the D2 and Di mesons is about 1.2, which falls within the wide uncertainty range of the 
current experimental value 1.3. The underprediction of the pionic decay widths of the L = 1 charm 
mesons is natural, as a substantial fraction of the total width is expected to be due to other decay 
modes, in particular tttt decay. The analogy with the corresponding decay modes of the K2{14:30) 
strange meson suggests that the decay modes D' -^ Dtttt may be responsible for a significant fraction 
of the observed decay widths. 

The results for the calculated pion decay widths of the excited D mesons obtained here are rather 
similar to those obtained in ref. |^, despite the different calculational framework and the different 
Hamiltonian model. The calculation in ref. |p| was restricted to the decays allowed by the lowest order 



selection rules suggested by heavy quark symmetry |10 . In the case of the Di meson the excluded 
S— wave transition was found to be rather important here, a result already hinted at in ref. |11[ . The 
conclusion reached here is that the chiral quark model does indeed provide a fair description of the pion 
decay widths of the orbitally excited D mesons, if they are treated as relativistic interacting two-quark 
systems. 

This paper falls into 5 sections. In section 2 the decay width for D* -^ Dtt is calculated. In section 
3 the corresponding decay widths for the orbitally excited Di and the D2 mesons, including the S-wave 
pion decay rate of the Di meson, are calculated. In section 4 the estimated pion decay widths of the 
hitherto undiscovered charm mesons with L = 1, J = 1 and L = 1, J = arc given. Section 5 contains 
a summarizing discussion. 



2 Pion decay widths of the V mesons 

2.1 Single quark approximation 

The main contribution to the decay widths of the D* mesons in the ground state band is due to the 
pion decays D* -^ Dtt. These transitions are intriguing in that since the mass difference Mw — Mo 
is very close to the pion mass, the available phase space is very small. Because of this closeness to the 
threshold for tt decay, and the nonzero mass splittings between the different charge states of the tt and 
D mesons, one of these decays - the decay of the D*^ to D^tt"^ - is in fact kinematically forbidden. The 
orbital wave functions of the constituents of the D and D* mesons differ very little from one another, 
which implies that the main pionic decay mechanism is P-wave pion decay. Because of the consequent 
threshold suppression and the small phase space, the total widths of the D* mesons are expected to 
be very small. The current empirical upper bound for the total width of the D*^ is 0.131 MeV and 
that for the D*'^ is 2.1 MeV [E2|. The former one of these upper bounds is already constraining for 
theoretical model calculations. 

The D and D* mesons are confined qc and qc systems, where q and q denote constituent m, d and 
■u, d quarks and c and c charm and anticharm quarks respectively. The pions only couple to the light 
flavor quarks, the simplest model for the coupling being the chiral coupling: 

g1 _ 

^ = iT^'^ql5lt,df,<j)„-T1pq. (1) 

Here g^ denotes the axial coupling constant of the light flavor constituent quarks and /tt is the pion 
decay constant (93 MeV). The value oi g\ should be somewhere between unity |1^ and g'^ = 0.87 |1J, 
p^ , [iql . Only the P-wave part of the coupling (p contributes to the pionic decays D* -^ Dtt, and 
hence, in the single quark approximation, the transition operator reduces to 

Tp = -J- u{p') j^j'' ■ k u{p) T^ 

g\ J E' + m J E + m / _ P^ - k^/A 



= -^27:V^^V^^l'"3(i?' + m)(P + m)j^'-^^'^- ^^^ 

Here m denotes the mass of the light constituent quark and k is the momentum of the emitted pion. 
The operator P is defined as (p' +p)/2, with p and p' being the initial and final momenta of the light 
quark. The energy factors E and E' that appear in eq. (g) are defined as y/p^ + m^ and v^j?'^ + m^ 
respectively. The pionic decay widths of the D* mesons may now be expressed as 

T{D*°^D\') = :^^ ('4-Xk'Ml (3) 

and 
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^ > 2A^Md.± \.fj "' ^ ' 

r{D*^-.D\^) =. T^W^i'fXk^Ml (5) 
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for the charged D* mesons. In the above expressions, A^o is the orbital part of the matrix element 
(00,00|CTg • fc |01, Im), where \LS,JM) denotes the state vector of the qQ system. To arrive at this 
expression, the following result for the sum over spins and integration over directions of k has been 
used: 

i /"dOfc^(lm|a, •fc|00)(00|CTg-fc|lm) = yfc^ (6) 

in 

Here ffq denotes the spin operator of the light constituent quark. The matrix element A^o may be 
expressed in unapproximated form as 
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TT 



Mo = - dr'r'uoir') / drruo{r) / dP P^ dz fssiP, z) 



"'0 "'0 J-1 



E' + m jE + m ( P^ -k^/A 



2E' \ 2E \ 3(S' + m)(S + m) 



Here uo{r) is the reduced radial wave function for the D and D* mesons. The explicit expressions for 
the energy factors appearing in eq. m) are 



E = y/m'^ + P2 _ pkz + fc2/4 and E' = y/m^ + P^ + Pkz + k^/A (8) 

respectively. The factor /bs {P, z) originates in the quasipotential reduction of the amplitudes defined 
for the Bethe-Salpeter equation, and is defined as 

M -\- m 



Here Ec = y/M'^ + P^ - Pkz + k'^/A and E'^ = a/M^ + P"^ + Pkz + k'^/A and M is the heavy quark 
mass. In the non-relativistic limit the expression (R) reduces to 

Mo-/ drulir)joi-y (10) 



Because of the small masses of the light constituent quarks the non-relativistic approximation (|1C 
turns out to be inadequate for a reliable description, and typically leads to overestimates of the exact 
result obtained from eq. (^ by at least a factor ^2. In order to evaluate the expression (|^), a 
model for the wavefunctions is required. In this work, the wavefunctions that were obtained by solving 
the Blankenbecler-Sugar equation in configuration space in rcf. |0| will be used. Those wavefunctions 
correspond to an interaction Hamiltonian which is a combination of relativistic one-gluon exchange and 
a scalar confining interaction. In this model the light quark mass was obtained as 450 MeV and that 
of the charm quark as 1580 MeV. The string tension in the scalar confining interaction was taken to be 
c = 1120 MeV/fm. The relativistic one-gluon exchange interaction employed in ref. [g also features a 
running color " couplant" , which was taken to be of the form 118] : 
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(11) 



In eq. (|ll|), Aq denotes the confinement scale, for which the value 280 MeV was used in ref. |^. For 
the gluon mass parameter mg, the value mg = 240 MeV was obtained. A comparison between the 
spectrum obtained in ref. B and the current empirical spectrum is shown in Fig. P. The resulting 
reduced radial wave functions for the IS"- and IP-states that were obtained with this model in ref. ||] 
are shown in Fig. 0. 



Empirical and calculated D meson spectra (State energies in MeV) 
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Figure 1: Experimental and calculated D meson states. The calculated states are solutions to the 
Blankenbecler-Sugar equation and correspond to the Hamiltonian used in ref. . The energies of the 
empirical states along with the pion momenta used in the calculations are listed in Table H. 



Decay 


Initial state mass [Mi) 


Final state mass (Af/) 


-K momentum (fc) 


i?*o ^ D°Tr° 


2007 


1865 


43.1 


D*± -^ D^TT° 


2010 


1869 


38.3 


D*^ -^ D°Tr^ 


2010 


1865 


39.6 


Dl -^ D*TT 


2422 


2009 


355 


D*^Dn 


2459 


1867 


505 


D^ ^ D*TT 


2459 


2009 


389 


Dl -^ D*TT 


2389 


2009 


326 


DI^Dtt 


2341 


1867 


408 



Table 1: Initial and final state D meson masses and emitted pion momenta in MeV used in the 
calculations. When different charge states are specified, the values are taken from [^, otherwise 
suitable averages are used instead. The D\ and Dq masses are taken from ref. Q. 
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Figure 2: The reduced wave functions for the S- and P-state qc (or qc) systems according to ref. p. 
As these wavefunctions correspond to the spin-averaged S- and P-states of the D meson, they do not 
incorporate the fine structure splittings of the various states in Fig. n^. Therefore, in order to avoid 
unnecessary errors, the masses and pion momenta listed in Table n are used in the calculations. The 
approximate wavefunctions corresponding to eq. (p5[) in Appendix A are shown for comparison. 



The numerical value of the unapproximated matrix element Ado, ^q. (M), as obtained using the wave 
function uo{r) shown in Fig. is M.o = 0.649. Although the value of the pion momentum k in the decays 
of the D* mesons is non-negligible (In the case D*"^ -^ D°n+ it is 39.6 MeV/c and in D*° -^ D°tt° 
43.1 MeV/c) the product of this momentum and the range of the wave function (~ 0.5 fm), is of the 
order ~ 0.1. Because of the smallness of this product the value of the non-relativistic approximation 



to the matrix element Mq as evaluated from eq. (10) is ~ 1.0. Thus one may conclude that even for 



the decays with the smallest pion momentum k, the non-relativistic approximation overestimates the 
calculated decay widths by about a factor 2.4. 

The calculated decay widths for D* -^ Dtt as obtained with g\ = 0.87 and g\ — 1 are given in 
Table y along with the pion momenta used and the present empirical upper bounds. The D meson 
masses used are those listed in Table ^ The calculated total pionic width of the D*^ is 0.082 MeV 
in the single quark approximation, a result, which is about 60% of the present empirical upper bound 
(0.131 MeV). In the single quark approximation the calculated width of the D*'^ is 0.036 MeV, which 
is far below the empirical upper bound of 1.3 MeV. In this case, the empirical upper bound is far too 
large for being theoretically constraining. The widths obtained here for the D* mesons are similar 
to those obtained in ref. |q], where the Gross equation framework was applied. In that reference, the 
value suggested for g\ was however only 0.75, and with that value the present calculated widths would 
be somewhat smaller. The difference may be attributed partly to the much lighter constituent quark 
masses used in ref. m and partly to the relativistic factors in eq. (R), which arise from the canonical 
boosts. Without those factors the calculated width of the D*^ would exceed the experimental upper 
bound in the single quark approximation. 

2.2 Axial exchange current contribution 

The single quark amplitude, eq. (0), represents a coupling of the pion to the axial current of the light 
constituent quark, which in the static approximation may be expressed as 

Aa = -g\SqTa. (12) 

The linear scalar confining interaction will contribute an exchange current that arises from coupling 
of the pion to a virtual negative energy quark, see Fig. 0. If the confining interaction is taken to be 
of the form Vc [r) — cr — b, as in ref. M , then the expression for the corresponding exchange current 
operator may, to lowest order in v/c, be written as p7[: 
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(13) 



This expression shows that the exchange current contribution is a relativistic correction of the same 
order in v/c as the relativistic corrections to the single-quark operator in cq. (0). The constant b was 
obtained as 320 MeV in ref. g. 

The expression (O) does not include the corrections from the canonical boost factors on the single 
quark spinors that are included in the single quark operator, eq. (0). Hence a more realistic evaluation 
requires that those factors are taken into account. For simplicity, the same spinor factors as for the 
single quark operator are used in this context, and the operator is also multiplied with the corresponding 
- numerically less important - spinor normalization factors for the scalar vertex on the heavy antiquark 
line. This implies going beyond the static local approximation for the confining potential. Accordingly 
the potential function cr is to be replaced with the general expression c\f' + r\/2. 



Moreover the factor l/m? in the axial exchange current operator (O) arises as the static approxima- 
tion to a combination of energy dependent factors appearing in the denominator for the intermediate 
negative energy state and of corresponding energy factors in the quark spinors. Consequently, the 
static approximation to these factors implies a very large overestimate of the axial exchange current 
contribution. In order to obtain a more realistic estimate for this contribution, the static approximation 
in the axial exchange current operator may be replaced as 



1 4 / 2m 

^ ^ 2m + E + E' [e + E' 



(14) 



The resulting matrix element of the exchange current operator (O) between the triplet and singlet 
states may then be expressed as 



■^S" = T^ I dr'r'uoir') I dr r uo{r) V^ ' "^ 




j'^dzfBs(P,z)jJr'^P^+'^ 



S{E' + m){E + m)^ 

(15) 

In the above expression, a symmetrized form \/{r"^ + r'^)/2 of the confining potential has been em- 
ployed for simplicity. This form leads to the correct limit when r = r' . Moreover, the term proportional 
to c/2r in eq. (|l^) has been neglected in the calculations because of its smallness. The numerical value 
of the exchange current matrix element A^q'^ as evaluated using eq. (|l^) is obtained as 0.040. In the 
static limit, i.e. without the replacement of eq. (|4|), that value increases by about a factor 3. The 
static limit is thus seen to represent a significant overestimate of this effect. In any case, the cor- 
rection is non- negligible in comparison with the contribution from the single quark operator, which 
is A^o = 0.649. Addition of the exchange current contribution increases the calculated pionic decay 
widths by about 15%, and brings them to the empirical upper limits in the case of the charged D* 
mesons, \i g\ — 1. These results are also shown in Table |[ 

Even though the total widths of the D*^ states have not yet been experimentally determined, their 
branching fractions for tt decay are known to a good degree of accuracy. The current experimental 
resuhs are (68.3±1.4)% for D*^ -^ D^tt^ and (30.6±2.5)% for L>*± -^ D^it° respectively |2). Using 
the results in Table 0, the ratio of these decay modes is obtained as 

being thus in excellent agreement both with the result obtained in ref. [B and the experimentally 
determined value 2.23 ± 0.19. 



Decay 


n momentum k 


SQA 


SQA + EXCH 


9', = 1 


Experiment 


n*± -> i^i^it 


38.3 


0.026 


0.029 


0.038 


<0.04 


D*± ^ D°7r^ 


39.6 


0.056 


0.064 


0.084 


<0.09 


D*° ^ i^o^o 


43.1 


0.036 


0.041 


0.054 


<1.3 



Table 2: The calculated pionic widths and experimental upper limits in MeV for the D* mesons, 
corresponding to g'^ = 0.87. The single quark approximation, with relativistic corrections is denoted 
SQA, and the result obtained with the exchange current contribution included is denoted SQA + 
EXCH. The net calculated widths are also shown for g\ = 1. 



3 Pion decay widths of the Di and D2 mesons 

3.1 Pion decay by the single quark axial current 

Only two of the four expected charm meson resonances with L = 1 have hitherto been discovered with 
certainty. These are the Di(2420) and the 152(2460) mesons, including all the different charge states. 
It is generally assumed that these are spin triplet states with J = 1 and J = 2 respectively. In the 
Heavy Quark Symmetry (HQS) framework, they are assumed to be states with light quark angular 
momentum jq = 3/2. For these D meson states the total widths have been experimentally determined, 
although with quite large uncertainty ranges, and some of their pionic decays have been "seen" [n2|. 
The basic pionic decay mode of these resonances is D— wave decay by pion coupling to the axial current 
operator, eq. (0), of the charm mesons. In the case of the Di meson, S*— wave pion decay through the 
axial charge operator also contributes significantly to the decay width pT[ . 

The axial current transition operator for pion decay to the ground state D mesons is given by 
eq. (0). If both the -Di(2420) and 1)2(2460) mesons are assumed to be mainly spin triplet states, the 
calculation of the pion decay widths of these states require the following spin sums: 
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for spin singlet final states, and 
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M=-J 
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for spin triplet final states. In the above expressions, 6*12 denotes the tensor operator 

Si2{k) = 3(Tq ■ kuQ-k- dq- dq. (19) 

Evaluation of the orbital matrix elements of these operators then leads to the following expressions for 
pion decay driven by the axial current operator: 

Ta{Di^Dtt) = 0, (20) 

for the Di meson, and 



TAD^^D-.) = _^(^^|.^m;, (21) 



i 2 Ed (g\\ ,3 ,,2 



VM^D,) = ^^~[f^) t-'Mi. (22) 

r^ni^D-.^ = ±||£l(£A)^=Mf. (23) 

for the D2 meson. In the above expressions, all final charge states have been included. A^i is an 
orbital matrix element defined in analogy with eq. ([7|) and may be expressed as 

-1 /"Oo /"Oc /"Oo />! 

Ml = - dr'r'uoir') / drrui{r) / dP P^ dz fssiP, z) 
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E' + m jE + m f P^-P/A \ k/A-Pz 



2E' M 2E \ 3{E' + m){E + m) J ^p2 + fc2/I6 - Pkz/2 
^o[r']IP^ + -, + — )n[r^P^ + -^- — j. (24) 

In eq. (p4), the factor /bs(^i 2) is defined as in eq. (O), and the quark energy factors E, E' are defined 
as for eq. (M). ui{r) is the reduced radial wave function for the D meson P-state 

111, JM) = J2 {nis\JM\ ^Yu{f) \ls) , (25) 

Is 

where |ls) denotes a spin triplet state with Sz = s. In the non-relativistic limit the matrix element (El) 
reduces to 

Mi=J druo{r)ui{r)jJ^\. (26) 

Note that eq. ( pO| ) predicts that the Di meson cannot decay to Dt:. This prediction only holds if 
the Di is a state with J = 1. As will be shown, it does not provide a means to discriminate between 
spin singlet and triplet states. However, the prediction that the D2 decays to both D* and D permits 

11 



the identification of the empirical D2(2460) meson as a spin triplet state with J — 2. The numerical 
values for the matrix element A^i, as obtained with the L = and L — 1 wave functions shown 
in Fig. for the pionic decay modes (pl[-p3[), are listed in Table & In all cases, the non-relativistic 
approximation overestimates the calculated matrix elements by ~ 50%, which is to be expected in 
view of the relatively small masses of the light constituent quarks. The calculated pionic decay widths 
corresponding to the matrix elements in Table S are given in Table a along with the current empirical 
values. 

3.2 Axial exchange current contribution 

The axial exchange current operator, eq. (|l3|), may also contribute significantly to the pionic decay 
widths of the L = 1 charm mesons. The orbital matrix element of the spin part of the exchange 
current operator between the L = 1 and L = states, when evaluated with the same spinor factors 
and approximations as in eq. (|lj), may be expressed as 



MT = TTT dr'r'uo{r') drrui{r)Vc 

127r Jo Jo 



dPP^ dzfBsiP,z) 

10 J -1 




„„ fc2 Pkz\. / 

^+16 + — FM^ 



k/A - Pz 



^2 k^ Pkz 

16 2 / ^P2 j^ fc2/i6 - Pkz/2 



2m 



E' \E + E' 



E' + m E + m 



lE'+M Er + M 



2E' 



2Ec 



2E' V 2E 

P^ - k'^/A 
iE'^ + M){Ec + M) 



P2 _ fc2/4 



3{E' + m){E + m) 



(27) 



The numerical values of these matrix elements are also listed in Table y. Their magnitude corresponds 
to ~ 10% of the matrix elements of the corresponding single quark operator, eq. (|2J), and evidently 
the exchange current contribution is numerically significant these cases as well. This contribution 
enhances the net calculated pionic decay widths of the L = 1 charm mesons and brings them closer to 
the empirical values. These results are shown in Table |5[ 



Decay 


SQA-NR 


SQA-REL 


EXCH 


Di -^ D*iT 


0.135 


0.093 


0.010 


Dl -^ Dn 


0.185 


0.126 


0.012 


D* -^ D*n 


0.147 


0.101 


0.011 



Table 3: Matrix elements A^i of the single quark axial current operator in the non-relativistic approx- 
imation, eq. (Pq), and without approximation, eq. (24), for pion decay of the Di and D2 mesons. The 
column "exchange" contains the matrix elements A^°^ corresponding to eq. (|27| ) of the spin part of the 
axial exchange current operator ([13|). 
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3.3 Pion decay by the axial charge 



The charge component of the axial vector coupling in eq. (|l|) also contributes to the decays of the 
-D- mesons with L = 1, and gives rise to both S- and I?- wave pion decay. The amplitude describing 
the coupling of the pion field to the axial charge component of the light constituent quark may then 
be obtained from eq. (nh as 



Ts 



9\ 



2m + E + E' 



P' +P 



(28) 



2/^ y/4:EE'{E + m){E' + m) 

where ui^ = ^Jk"^ + m^ is the energy of the emitted pion. In the evaluation of the matrix element of 
Ts between Z3-meson states with L = 1 and the ground state charm mesons, the operator (p ' + p ) /2 
is treated as the differential operator i (V ' — V)/2. Application of the operators on the initial and final 
orbital wave functions leads to spin matrix elements of the form 



{S'M\ J2 {nis\JM) a] lis) = -5js'5m'm [VsJs'o + V25s'i 



(29) 



from which it follows that iS-wave pion decay can only contribute to the Di — > _D*7r decay widths. The 
resulting contribution to the -Di(2420) — > _D*7r decay width from S'-wave pion decay may be written in 
the form 



as 



Tc{Di 



D*7r) 



1 Ed> (g\ 



327: Md, \U 



■AV I^Y kM' 



IS- 



(30) 



In eq. (pO[), the matrix element Alis is defined, in analogy with the matrix elements in eqs. (R) and (24) 



M 



IS 



-1 /"Oo /"Oo /"Oo /•! 

- / dr'r' / drr dP P^ / dz fssiP, z) 

"' Jo Jo Jo J -I 



m{2m + E + E') 
^AEE'{E + m){E' + m) 

jo 



Mo(r')Mi(r) - u^{r')u'^{r) 



■uo(r')wi(r) 




p2 ^ 

16 2 



(31) 



In the above expression, the factors u'{r) denote the derivatives of the reduced radial wavefunctions 
displayed in Fig. y. In the non-relativistic limit, the matrix element ( pl[) simplifies to 



Mis = 



dr 



Uo(r)wi(r) ~ UQ{r)u'-^{r) 



uo{r)ui{r) 



kr 



(32) 



The numerical values of these matrix elements are given in Table 0. The values of the matrix elements 
Alls in Tabic are very large and lead to unrealistically large contributions from S'-wave pion decay 
to the widths of the Di mesons. Insertion in eq. ( pO[ ) would give for Tc [Di -^ D*t:) the value ~ 60 
MeV, which exceeds the empirically determined width by about a factor 4. Reduction of this value to 
a realistic level is however brought about by the two-body mechanism described in Section [3.4| that is 
implied by the linear confining interaction, if its spinor structure is scalar. 
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The amplitude (|2^) also gives rise to -D-wave pion decay, which contributes to the decay widths of 
both the Di and Z?2 mesons. The explicit expressions for these contributions from D-wave pion decay 
may be obtained as 



Td {Di ^ D*7r) 



1 Ed' 



SttMd, \U 



for the Di meson, and 



^d{d; 

Td {d; - 



■^ Dn) = 
D*tt) = 



9\ 



\ m J 



kM^ 



ID 



(33) 



i 2 Ed 
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8tt 5 Afoj 
3 3 Ep, fg\ 
8tt5Md* \fn 



\ m J 



kM 



IDi 



\ m / 



ID- 



(34) 
(35) 



for the D| meson. In the above expressions, the matrix element AIid is again defined analogously to 
eq. (pi|), and may be expressed as 



M 



ID 



-ir'^-'i?"r''''i'^'-^^'->(ii^T0s 



m{2m + E + E') 
y/'iEE'{E + m){E' + 



Uo(r')Mi(r) - UQ{r')u'^{r) 



, , , ^ fc2 Pkz \ / / ^ /c2 pk-^ 



UQ{r')ui{r) 



(36) 



In the non-rclativistic limit, this expression reduces to 



M 



ID 



dr 



Uo(r)Mi(r) - uo{r)u'-^{r) 



ua{r)ui{r) 



kr 
J2ly 



(37) 



The numerical values of these matrix elements for Z?— wave decay as obtained with the wavefunctions 
of ref. ||6|] turn out to be exceedingly tiny in comparison with those from S-wave decay. The reason 
for this is immediately apparent if one considers the wave function combination u'^ui — uou[ + uoUi/r 
that appears in the integrand in eq. (|3^). If the approximate wave functions plotted in Fig. Q and 
given in Appendix A are used, that wave function combination vanishes exactly. As these approximate 



wavefunctions are not perfect and as the nonlocal structure of the integrand in eq. (36) has to be taken 



into account, this cancellation is not absolute if the unapproximated expressions and wave functions 
are used. The numerical values of eqs. (|6|) and ( p7| ) turn out to be of the order 1 MeV and thus 
completely negligible as compared to the matrix elements of S*— wave decay listed in Table H. These 
I?-wave decay amplitudes that arise from the axial charge operator also obtain a contribution from the 



two-body mechanism outlined in Section 3.4 
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3.4 Two-quark contributions to the axial charge operator 



Both the confining and one-gluon exchange (OGE) interactions contribute two-body terms to the axial 
charge density operator. If the linear confining interaction has the form Vc{r) S, where Vc is the central 
confining potential, and S is the scalar (Fermi) invariant for the two-quark system, it gives rise to an 



axial exchange charge density operator, which to lowest order in v/c may be expressed as |19 






Pr„ 



The corresponding expression for the vector coupled one-gluon exchange interaction is 



^OGE 



mM 



Vg{r) a ■ P-,Ta, 



(38) 



(39) 



where the momentum operator Pg is defined as Pg = (pe '+Pg)/2. Here Pc and pg ' denote the momentum 
operators of the heavy quark and Vg (r) is the main spin- independent term in the one-gluon exchange 
potential. If the color couplant as is taken to be constant, then Vg{r) may be expressed as Vg{r) = 
—Aas/Sr in the static limit. 




K 




K 



Figure 3: Two-quark contributions to the pion production amplitude associated with the axial charge 
operator. The diagrams shown correspond to both time orderings of the two-quark contribution from 
the scalar confining interaction. Two similar diagrams arise from the one-gluon exchange interaction, 
in which case the scalar vertices are to be appropriately replaced. 



In the static limit, the two-body contribution that arises from the confining and one-gluon exchange 
interactions to the amplitude for S'-wave pion decay may be expressed in the form 



t: 



(2; 



1 / 4: m as 

= \ cr — b -\ 

TO V 3 M r 



Ts. 



(40) 
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Here, Tg denotes the single quark amplitude for S'— wave pion decay, eq. (|2^). This result may be 
derived as a pair current operator, or more simply by making the scalar shift in the light constituent 
quark mass m -^ m + Vc{r) in the case of the scalar confining interaction. This two-body operator 
arises in the non-relativistic reduction of the general pion decay amplitude for the qQ system as a 
pair (or seagull) term, see Fig. ^. A corresponding two-body term appears in the amplitude for pion 
production in nucleon-nucleon collisions, the difference being that K is in that case replaced by the 
effective scalar component of the nucleon-nucleon interaction. 

In eq. (40) the factor 1/m represents the static limit of the propagator of the intermediate negative 
energy quark. As the static limit is known, from the treatment of the axial exchange current, to give 
large overestimates in the case of light constituent quarks, the static propagator will here be replaced 
by the symmetrized form 4/ {2m -I- £^ -I- E'). The matrix element of the two-quark contribution to the 
S*— wave decays through the axial charge operator may then be obtained by modifying the single quark 
matrix element for S'-wave pion decay, eq. (|3l|), according to eq. ([40|), giving 



M 



2S 



1 



dPP" 



dzfBs{P,z) 



Am 



^/AEE'{E + m){E' +m) 
jo I r 



u'oir')ui{r)-Uo{r')u[{r)-2 



Uo{r')ui{r) 




(41) 



In the above expression, the symbol ~ indicates that |r+r '|/2 has been approximated by the expression 
■\/(r'2 -I- r2)/2 as in eq. ([l5|). In the nonrelativistic limit, eq. ( [4l| ) becomes 



X™s"f 



- / drVcir) 



Uo{r)ui{r) - UQ{r)u[{r) 



uo{r)ui{r) 



kr 
Joly 



(42) 



The gluon exchange contribution, eq. (Ba), to the axial exchange charge operator gives a contribution 
of smaller magnitude than the confining interaction because of the heavy antiquark mass involved. Since 
the heavy quark constituent mass, 1580 McV, is not very large compared to the light quark mass of 450 
MeV, the onc-gluon exchange contribution turns out to be significant as well, see Table ^ and cannot 
be neglected. It contributes the following orbital matrix element for pion decay of the Di meson: 



A^?s°^ 
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7^ JO JO 

Am^ 



drr Va 




dPP' / dzfBs{P,z) 
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Uo(r )ui[r) - uoir )ui(r) - 2 
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Jo r 




(43) 



The effective one-gluon exchange potential Vg{r) that is appropriate for the D meson systems and which 
takes into account the screened running color coupling as has been calculated in ref. [g] . This potential 
function may be well parametrized as Vg{r) — —A a,rcta,n{Br)/r, with A=1.1899 and B=3.39768/fm. 
Using this form, the numerical value A^^g^^= 132 MeV is obtained. In the static non-relativistic limit, 
the one-gluon exchange contribution reduces to 



M 



OGE 

2S 
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M 



drVgir) 



Uo(r)Mi(r) ~uo{r)u[{r) 



uo{r)ui{r) 



30 



(44) 



The two-quark mechanism shown in Fig. contributes to the D-wa.ve pion decay amplitude from 
the axial charge operator as well. As in the case of S'-wave pion decay, this effect may be taken into 
account by adding to the one-body matrix element A^id the matrix element that is associated with 
the exchange charge contribution. For the confining interaction, this matrix element may be obtained 
similarly to eq. (BTl), giving 



>is™f 
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(45) 



A similar expression may readily be constructed for the Z?— wave contribution from the one-gluon 
exchange interaction by the techniques outlined in this section. The numerical value of the matrix 
element ( p5| ) is of the same magnitude (1 MeV) as that of eq. (^6|) and thus completely insignificant. 
As the one-gluon exchange contribution is smaller by a factor m/M , it will not be considered in this 
paper. In the nonrelativistic limit, the expression (MS) reduces to 



A^conf 



1 



drVc{r) 



u'^{r)ui{r) -uo{r)u[{r) 



UQ{r)ui{r) 



n 



(46) 



The numerical values of the matrix elements for S'-wave pion decay of the Di meson are given in 
Table 0. The matrix element of the two-quark contribution A^2s"^' ^"^^ (E3)' ^^ ^^ ^^^ same magnitude as 
the single quark contribution, eq. (pl|), and has opposite sign. The corresponding matrix element from 
the one-gluon exchange interaction, eq. (Eq), represents a somewhat smaller contribution. Addition of 
these three matrix elements reduces the net matrix element for S'-wave pion decay of the Di state by 
about a factor 3 as shown in Table 0. With this smaller value the net contribution of S-wave pion 
decay to {Di -^ D*tt) becomes quite small and comparable to that from the axial current operator, as 
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indicated in Table H. It should, however, be noted that the results are very sensitive to the exact form 
and composition of the interaction Hamiltonian used in the calculations. 

In the case of the D2 mesons the calculated width falls below the current empirical range by about 
40% if g\ = 0.87. With g\ = \ the calculated values are only slightly below the empirical values. 
These results leave room for about 5-10 MeV for the contribution from tttt decay. This is what would 
be expected on the basis of analogy with the decay pattern of the strange iir|( 1430) meson, which 
should have a structure similar to that of the I?2(2460), once the charm quark is replaced by a strange 
quark. If the nonrelativistic values for the orbital matrix elements M.\ given in Table [^ are used in 
the calculation, the calculated width for the D\ meson would be well above the empirical range. The 
ratio of the calculated pion decay widths of the D\ meson and the Di meson is 1.2. This result is 
compatible with the current empirical value ^ 1.3. 



Decay 


A^is 


X™"f 


A^OGE 


Total 


Di ^ D*t:, NR 
Di -^ D*TT, REL 


-1194 
-592 


+450 

+278 


+411 
+ 133 


-333 

-181 



Table 4: Matrix elements in McV of the single quark and two-quark operators for the S*— wave axial 
charge contribution to the decay Di —^ Dn. The resulting net contribution to this decay mode is also 
given. The labels NR and REL indicate that the non-relativistic and relativistic expressions have been 
used respectively. 



Decay 


Current 


Charge 


Total 


9'a-I 


Experiment 


Di -> D*t: 


4.2 (3.4) MeV 


6.1 MeV 


10.3 MeV 


13.6 MeV 


18.9l|:| MeV 


D^ -> Dtt 


8.1 (6.7) MeV 


- 


8.1 MeV 


10.6 MeV 


7 


DI -> D*t: 


3.9 (3.1) MeV 


- 


3.9 MeV 


5.1 MeV 


7 


D^ -^ Dn + D*TT 


11.9 (9.9) MeV 


- 


11.9 McV 


15.7 MeV 


25l? MeV 



Table 5: Calculated and empirical pion decay widths of the Di and D2 mesons driven by the axial 
current and charge operators respectively, for g\ = 0.87. The empirical values are total widths p2[ , 
which should mainly be due to pion decay to the ground state. The numbers in parentheses are the 
decay widths obtained without the axial exchange current contribution. The calculated values are also 
shown for g\^ \. 

In the case of the Di mesons the total width also obtains a significant contribution from S'-wave 
pion decay. Here it has been assumed that the empirical widths are mainly due to pion decay to the 
ground state D and D* mesons. In addition, even though the present empirical data on the L — 1 
D mesons is severely limited, the ratio of Dt: to D*tt decay of the D2 meson has been measured [ [I2| . 
From the results in Table |5|, we obtain 
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r [Dl -^ Dtt) 

r [Di -> d*tt) ~ ^•^' ^^^^ 

which is in agreement with the experimentally determined value 2.3 ±0.6 for the neutral Dg meson |1^ 
(The measured value for the charged D2 meson has much larger statistical uncertainties). Note that 
in Table |5|, the quoted experimental values are for the neutral Di and the charged D2 mesons. The 
corresponding value for the neutral D2 meson is 23 ± 5 MeV |1^ . The available data for the charged 
Dl meson is considerably poorer, since it has only recently been discovered. 



4 Pion decay widths of the D mesons with L = 1, J = 0, 1 

Two D meson resonances with L = 1 remain to be found experimentally. In the LS basis these are 
the spin singlet state Dl with J — 1 and the spin triplet state Dq with J = 0. In this context, it is 
worth noting that although the 'star' in the labeling notation D* is usually reserved for states with 
J^ = 0+, 1~, 2+, • • •, it has become conventional to label the spin singlet state Dl to distinguish it from 
the spin triplet state Di. In the Heavy Quark Symmetry (HQS) framework, the D^ and Dq mesons 
correspond to states with light quark angular momentum jq = 1/2. Both of these states may pion 
decay by D—wave decay through the axial current operator as well as by S*— wave decay through the 
axial charge operator. Consider first the (mainly) singlet Dl state. Because of the spin dependence 
of the pion-quark coupling, it follows that this state can only pion decay to the triplet D* meson. To 
derive the expression for the width of this decay mode, the following additional spin sums are required: 

Ss ^ i; Yl (10,lA/|ag-fc|00,00)(00,00|CT, •fc|10,lM) 

1 ^ 1 ^ 

= - ^ (10,lA/|--5i2(fc)|10,lAf), (48) 

M=-l 

for spin singlet final states, and 

ill 
^* = ^ Yl Yl (10,lAf|(7q-fc|01,lTO)(01,lm|aq-fc|10,lM) 

A/=-lm=-l 
1 1 1 ^ 

= - 5] (10,lAf|fc2 + -5i2(fc)|10,lM), (49) 

for spin triplet final states. From these expressions, it follows that the decay mode Dl -^ Dtt is 
forbidden, since the tensor operator 6*12 has a vanishing matrix element between spin singlet states. 
The resulting expressions for the pionic decay width of the DJ state driven by the axial current operator 
is then 
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The corresponding expression for the pion decay width of the spin triplet Dq state may, using the spin 
sums given in eqs. (|l^) and (|l8|), be expressed as 



r^ (D* ^ Dn) 



3 Ed 



9l 



87rA/n. V/- 



k^Ml 



(51) 



The axial charge operator also contributes to these decay modes. The expressions for these 5-wave 
contributions may be obtained as 



for the Dl meson, and 



Tc {Dl -> D*7t) 



Tc {D* ^ Dn) = 



6Att Md 



'"■ 'f ) (^) ^^;.' 



3 E 
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647rMn. \U 
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(52) 



(53) 



for the Dq meson. The matrix elements required for the evaluation of these decay widths are given in 
Table y, and the pion momenta and masses used can be found in Table |l|. The resulting calculated 
pion decay widths of the D\ and Dq mesons are given in Table 0. It is evident that the S-wave 
contributions are dominant in these cases, although the D-wave decays also contribute significantly. 
In the HQS framework, only S'-wave pion decay is allowed to contribute H. 



Decay 


Ml 


MT 


Alls 


Mir' 


XopE 


A^™T 


Dl -> D*t: 
D*Q -^ Dtt 


0.086 
0.106 


0.010 
0.011 


-598 MeV 
-578 MeV 


-1-284 MeV 
-f 267 MeV 


+135 MeV 
-hl30 MeV 


-179 MeV 
-181 MeV 



Table 6: Matrix elements of the axial current and charge operators for the decays of the D* and Dq 
mesons. The matrix elements Mj^"^ represent the resulting net contribution to S'-wave pion decay 
from the axial charge operator, eq. (|2^), when the two-body matrix elements from Section B.4 
added to the single quark contribution. 



are 



Decay 


Current 


Charge 


Total 


9'a = 1 


Dl -> D*Tr 


2.8 (2.3) MeV 


7.2 MeV 


10.0 MeV 


13.2 MeV 


D*Q ^ Dtt 


7.9 (6.5) MeV 


13.0 MeV 


20.9 MeV 


27.7 MeV 



Table 7: Predicted pion decay widths of the Dl and Dq mesons driven by the axial current and charge 
operators respectively, for 5^ = 0.87. The empirical values are total widths |12| , which should mainly 
be due to pion decay to the ground state. The numbers in parentheses are the decay widths obtained 
without the axial exchange current contribution. The calculated values are also shown for g^ = 1. 
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The pion widths of the D^ and Dq states have been predicted to be much larger in the previous 
hterature ^, ^, although in ref. ^ a considerable reduction of these widths was already hinted at. 
One reason for the smallness of the values obtained here is the more complete degree of suppression 
of S— wave pion decay modes achieved here. This suppression was already shown to be necessary in 
order to avoid a large overprediction of the width of the Di meson. Since eqs. ( pO| ) and ( p^ ) suggest 
that 

Tc {Dl ^ D*7r) = ^ Tc pi ^ D*7r) , (54) 

if the pion momenta k are the same, it immediately follows that the S'-wave widths of these two states 
are of the same order of magnitude. Consequently, as the triplet Di state does not have any large 
S'-wave contribution to its pionic width, it follows that the singlet Z?* state should be rather narrow 
as well. The calculated pion decay widths of the D* and Dq mesons are here further reduced by the 
significant spin-orbit splittings in the P-shell, see Fig. |l|. The spectrum obtained in ref. |^ predicts 
the traditional " hydrogen- like" ordering of the P-wave states, which is supported by recent NRQCD 
lattice studies M. Since the decay rates are very sensitive to the amount of phase space available, this 
effect turns out to be quite significant. Consequently, the D2 meson has considerably more phase space 
available than the other P-shell states. 

In this context, it is instructive to note that the pion decay rates of the Kq{1430) and iir2 (1430), 
which are analogous to the Dq and D2 states, have been experimentally determined jlj]. The ratio 
of the widths of the Ktt decays of these strange mesons is particularly interesting, since the pion 
momenta in these decays are exactly equal. The ratio of the width of Kq — > Ktt to that of K2 ~> Ktt 
is ~ 5.7. Using eqs. (O) and (p3) with the pion momentum and matrix elements of the D2 meson, the 
corresponding ratio for the D mesons comes to ~ 5.4, which supports the conclusion that the pionic 
width of the Dq state should not be much larger than ~ 30 MeV, as indicated in Table M. 

It should however be noted that as a consequence of the cancellation between the exchange charge 
operator contribution and the contribution from the single quark axial charge operator, the S'-wave 
decay widths are very sensitive to the strength, form and composition of the interaction Hamiltonian 
of the qQ system. Empirical determination of the decay widths of the Dq and Dl mesons should as a 
consequence be able to provide useful information on the forms of these interactions. 



5 Discussion 

The present calculation of the pion decay widths of the excited charm mesons indicates that the chiral 
quark model description is compatible with the extant experimental data. The underprediction of the 
total decay widths of the L — 1 charm mesons is to be expected, as the pion decay pattern of the 
strange meson L — 1 states as e.g. the -ftrj(1430), suggests that tttt decay should give a contribution 
to the decay width of the D2 meson that is about a third as large as that of single pion decay. 

This point illustrates the difference between the present relativistic Hamiltonian method and the 
heavy quark effective field theory (HQET) method, in which at the present stage the effective coupling 
strength of pions to light constituent quarks has to be determined by the empirical width of one of the 
orbitally excited D mesons under the assumption that the width is solely due to single pion decay [gO| . 
Given that restriction, the present results for the decay widths of the Di and Dj mesons are overall 
rather similar to those in ref. pOl, although some details as e.g. the ratio between the decay widths to 
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Dtt and D*Tr of the D2, eq. (p7|), obtained there is 3.5 and thus somewhat above the empirical value 
2.3 ±0.6. 

The present calculation, which was carried out within the framework of the Blankenbecler-Sugar 
reduction of the Bethe-Salpeter equation, affirms the observation of ref. S], that a relativistic treatment 
of the qQ system is required, as a realistic description of the pionic decay widths is not attainable 
with a non-relativistic description. The framework of the Blankenbecler-Sugar equation brings the 
advantage of formal similarity to the conventional non-relativistic quantum mechanical treatment. In 
this approach the exchange currents that are associated with the scalar confining and vector one- 
gluon exchange interactions play an important role in suppressing the otherwise unrealistically large 
amplitude for S*— wave pion decay of the Di meson. This also has the effect of drastically reducing the 
calculated widths of the singlet D^ and triplet Dq states, which are predicted to be extremely broad 
in the HQS framework. In that framework, the allowed pion decay modes are determined by the total 
angular momentum of the light constituent quark, which is appropriate if the spins of the light and 
heavy quarks are decoupled. This is the case if the heavy quark mass is very large compared to that 
of the light constituent quark. However, the empirical D meson spectrum shown in Fig. |l| indicates 
that the spins of the quarks are strongly coupled, as the D* — D splitting is of the order ~ 130 MeV. 
In view of this, the present treatment is more general since it is not restricted to the case when the 
heavy quark mass is infinite. The present work incorporates all decay modes that contribute in that 
limit m. An analogous treatment within the framework of HQS requires an l/Mg expansion along 
with the introduction of new parameters, which have to be fitted to experimental data. 

The main goal of studying the pion decays of the charm mesons is that of determining the numerical 
value of the axial coupling g^ of the light constituent quarks. As long as the absolute values of the 
widths of the D* mesons are not known, this goal cannot be conclusively attained. Nevertheless both 
the present study, as well as the earlier investigation in ref. ||^, suggests that the standard values of 
g^, being somewhat below 1 appear to apply well in the case of the pion decays of the charm mesons. 
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Appendix A 

Because of the smooth radial dependence of the scalar confining and the one-gluon exchange interactions 
in the qQ system, the reduced wave functions of the low-lying states of the qQ system may be well 
approximated by the following simple expression: 

ui{T)^Nir'e-'''-"'\ (55) 

Here rrir is the reduced mass of the qQ system and c is the confining string tension. The coefficient d 
is defined as 

2 

d = -^/2mrC, (56) 

and Ni is a normalization factor chosen so that 

/•OO 

/ dru/(r)2 = l. (57) 

Jo 

These approximate expressions have the correct behavior at both small and large values of the quark 



separation r. For / = 0, 1 the expressions for the normalization constants may, with the aid of eq. (56) 
be obtained as 



No = -y/inirC, 

The approximate reduced wave functions uo{r) and ui(r) are compared to the corresponding numerical 
solutions to the Blankenbecler-Sugar equation in Fig. 
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